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'nI" ' ABSTRACT 

o ; 

' At the LHC, superpartners with the masses hghter than a few TeV may be found, 

and the masses of the supersymmetric (SUSY) particles and their interactions 



(N 



' will be studied. The information will be the base to consider the SUSY breaking 



(N 

I 



mechanism and cosmology related to the SUSY dark matter. We review the 
recent studies that aims for the "precise SUSY study" at the LHC. 



o 

1. Introduction 

> . 

. The LHC experiment, scheduled to start from 2007, is the best place to study the 

supersymmetry in coming 10 years. The LHC is a pp collider at ^/s = 14 TeV. It is able 
to find q and g in the minimal supersymmetric standard model(MSSM) up to 2-3 TeV. 
Not only for that, it is also the place to do the precise measurement of the nature of the 
O [ SUSY particles. If the sparticle masses are as light as 1 TeV, the production cross section 
is as large as 3 pb, therefore 0(10000) events /year will be produced for the detailed 
study. It is our responsibility to extract physics information, as much as possible, from 
the LHC experiment for the future developments of the experimental and theoretical 
^ ! particles physics. 

Different SUSY breaking models leave special imprints to the mass spectrums, there- 
^ ; fore measurement of the masses of the SUSY particles is a important target for the collider 
^ \ physics. After finding superpartners, we can select pure samples of the SUSY events free 
from standard model (SM) backgrounds. By looking into the distributions of the sparticle 
decay products, we can determine the mass of the sparticles, which may provide us the 
understanding of the SUSY breaking mechanism, namely, physics at the very high scale. 

The measurements also have impacts to the other physics. The mass spectrum deter- 
mination constrains the SUSY contributions to the fiavor violating processes, such as b 
rare decays, fi — > ej, and muon anomalous magnetic moment. The thermal relic density 
of the SUSY dark matter also can be calculated if the MSSM parameters are known. 

One may even study some of the SUSY relations. In the MSSM, the chiral structure of 
the coupling of the superpartners are restricted. For example the bino and wino couplings 
are the (s)fermion chirality conserved form fL(R) fL{R)B or /lIlW. The sparticle left-right 
mixing are induced by the F term as {m'^)j^ oc m^/i tan/?. We will discuss how to prove 
such relations at the LHC for explicit examples. 

Before going into the detail, let us first summarize basic features of the collider sig- 
natures of the SUSY particles at the LHC. In the models where SUSY is broken at high 
energy scale, the squark and gluino masses are much higher than that of the weak inter- 
acting SUSY particles. Therefore, at least one jet originated from a squark or a gluino 



cascade decay has large px- The decay products include weakly interacting SUSY parti- 
cles, which further decay into the leptons or jets. Especially, at the end of each sparticle 
cascade decay, the lightest SUSY particle (LSP) appears, which may be stable due to 
the R parity conservation. The LSP cannot be detected in the detector, therefore SUSY 
events have large missing momentum. By applying the cuts on the missing and the 
total transverse energy M^s, the SM backgrounds can be reduced to the negligible level. 

One of the relatively clean channels is the 2 body cascades or 3 body decays of the 
heavier neutralinos which have been obtained lots of attentions in the past physics 
performance studies [T], 

X° - (//) - llx'v (1) 

The cascade decay produce opposite sign same flavor leptons(OSSF). The accidental lep- 
tons can be estimated from the e/i distribution (opposite sign opposite flavor leptons — 
OSOF), leaving very clean samples containing 

Another important SUSY channels are those containing b quarks. In many SUSY 
breaking models, b and t could be much lighter than other squarks because of the negative 
running of the squark mass by the Yukawa coupling, in addtion to the left-right mixing. 
The decays g ^ b,i are to open in the wide parameter space. The LHC detectors have a 
good tagging efficiency for b jets ~ 0.5, and the cascade decays, 

maybe studied at the LHC by tagging two b jets. In the next section, we will discuss 
these processes as well. 

2. Sparticle Mass Measurements at the LHC 

If the LSP is stable, the SUSY events have large missing transverse momentum. This 
is important signature to discriminate it from the SM backgrounds. On this other hand, 
one has to determine the parent SUSY particle masses from the momentum of the visible 
particles only. This is easy for the case of e~^e~ collider experiments. The energy of 
collisions are known, and the decay kinematics can be solved by the constraint. But it is 
not simple task at the LHC, where a proton is a composite particle and collision energies 
of the partons are not fixed. 

The techniques to do the jobs at the LHC have been developed in the past 10 years. 
We quickly review the previous works and also explain new ideas. 

2.1. The case with the long lived NLSP 

In the models with the gravitino LSP, the NLSP is long lived and could be either 
charged or neutral. The NLSP candidates are sleptons (especially the lighter stau fi), the 



lightest neutralino Xi? or sneutrinos u. The NLSP could decay into its SM partner and a 
gravitino 

Xi 7^3/2- (3) 

One can extract the mass of SUSY particles by measuring the endpoints of the invariant 
mass distributions of the selected particles involving the photon or lepton. (The endpoint 
method is discussed in the next section). For example, from the endpoint of the m^i, 
m^ii and mu .. of the cascade decay X2 ~^ ^ Xi^^ ^ ^3/2^^7, one can calculate Xi(2) 
and I mass model independently with high precision of 0(1) GeV. The structure of the 
messenger sector will be determined very precisely from the mass measurements. 

The life time of the NLSP depends on the scale of the hidden sector SUSY breaking 
F which is also related to the gravitino mass 1123/2 = F/{\/3Mpi). We may learn the 
structure of the hidden sector from the measurement of the life time. For example, the 
life time is given for the decay / I1P3/2 as 

/100GeV\^ f Y 1, 

rNLSP= T7^^ 3x10-1^ 4 

V mNLSP / \lOOTeV J 

The gravitino could be a dark matter. It can be produced either thermally or from 
the NLSP decays. However if the life time of the NLSP is longer than 1 sec, the decay 
products may affect the BBN and effect depends on the life time strongly. The cosmology 
related to the gravitino is discussed in this workshop in detail |H3] . 

The signature of the charged NLSP(CNLSP) would be very spectacular. It is detected 
as a highly ionizing particles as it is non-relativistic. The time of the flight (TOF) infor- 
mation at the muon systems and momentum measurement at the inner detector allows 
precise measurement of the NLSP mass[3]. All neutralino masses may be determined 
within the error of 0(1) GeV by combining the NLSP momentum with that of the other 
visible particles. 

The supergravity interaction may also be studied in detail for this case. The charged 
NLSP may travel through the detector with finite vertex until it decays finally. The 
life time can be calculated from the distribution of the flight distance and the decay 
time of the NLSP. It should be possible to do such measurements as far as they decays 
inside the detector, although there are no systematic experimental study on the expected 
precision of the lifetime at the LHC. Even if ctnlsp is longer than the detector scale, the 
stable charged particle may be stopped at the massive stopper placed nearby the main 
detectorjl). Recent study shows that one can determine the life time as long as O(IOO) 
years at LHC, if a massive stopper with mass around 1 kton can be placed near the LHC 
detectors |6l7j . 

For the neutral long lived NLSP, the NLSP momentum cannot be measured directly. 
If the life time is short enough, it decays into a photon and a gravitino in the detector, and 
the photon will be detected at the ECAL. The reconstruction of the NLSP momentum, 
a decay position and a decay time of the NLSP is possible event by event by using the 



time information and photon momentum measurement at ECALjS]. This is because non- 
relativistic neutrahnos fly in the detector then decay, therefore the arrival time depends 
on both the NLSP momentum and decay time. The resolution of the decay time and 
position depends on the time resolution of the ECAL. The ECAL time resolution is 
about O(O.l) nsec in a beam test, which corresponds to 0(3cm) resolution on the NLSP 
vertex. Assuming the time resolution, the mass error obtained from the reconstructed 
NLSP momentum is estimated as low as 0(1) GeV at a model point jH], for the process 

I ^ Xll ^ 4^s/2ll. (5) 

In the limit where the NLSP life time is too long to decay in the detector, the signature 
becomes very close to the case where the LSP is the lightest neutralino. This will be 
discussed in the next subsections. 



2.2. The case with the lightest neutralino LSP at the LHC 

When the lightest neutralino is stable, it is not possible to measure the momentum 
directly. The study of the mass spectrum is still possible by looking into the decay 
distribution of the visible particles from sparticle decays. 

The past physics studies are mostly done with the endpoint method. Namely,one 
measures endpoints of the invariant mass distribution of the selected sample which are 
dominated by the events from a decay cascades. The positions of the endpoints depend 
on the sparticle masses. For example, the endpoint of the m{l^l~) distribution for the 
cascade decay given in Eq. ^ is 

mii{3 body) = m^o — m^o, 



mii{2 body) 



\ 



(m|o - m?)(m? - m|o) 



(6) 



respectively. One can then select a jet from a squark decay q qy^ by taking one of the 
two highest pt jets. The endpoint of the invariant mass of the jet and leptons provide 
the constraints among rriq and the lighter SUSY particle masses involved in the cascade 
decay. 

For the point SPSla which is defined in the mSUGRA scenario by the parameters 
tan/3 = 10, m = 100 GeV, M = 250 GeV, Aq = -100 GeV, and /i > [9;, the error of 
sparticle masses at the LHC for the integrated luminosity / Cdt = 300 fb~^ are estimated 
in [TUj, and the results are summarized in Tabled 

The cross section at SPSla is large (cr ~ 58pb) and O(IO^) sparticle will be produced 
at the high luminosity option of the LHC at this point. For weakly interacting particles, 
the errors are dominated by the uncertainty of the absolute mass scale, while the mass 
differences of the sparticles are more strongly constrained. This is general feature of 
the determination of the sparticle masses when the LSP momentum cannot be measured 
directly. The mass of the lightest neutralino may be measured more precisely at the 
international linear collider (ILC), then error of the sparticle masses would be reduced 
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Table 1: Mass of the sparticles at SPSla, and expected sensitivity at the LHC, and that with the expected 
input from the ILC. 



significantly as can be see the 4th column of the Tabled Finally the errors for the squark 
and gluino masses are dominated by the 1% jet scaling error. 

It is notable that the LHC can access to the mass of the the heaviest neutralino X4 
through its decay into sleptons. For this model m^o = 378 GeV, which cannot be accessible 
at ^/s = 500 GeV e~^e~ colliders. The measured mass difference m^o — m^o directly 
constrains fi parameter. The errors for the MSSM parameters would be significantly 
reduced when the measurements at the LHC and the ILC can be combined jTJ. 

We also can see the implication to the cosmology. The mass of the slepton and the 
lightest neutralino can be determined within 0(5)% level. This means one can calculate 
the neutralino pair annihilation rates into the first and second generation leptons at the 
time of decoupling with in the error of 10%. The pair annihilation rates into f is also 
important. The f mass also can be measured at the LHC, but more study is needed to 
obtain the error, as the channel involving r leptons suffers backgrounds from QCD jets. 



2.3. Mass relation method 



Although the mass determination through the endpoint method is successful, there are 
problems which may limit the application. The problems may be summarize as follows; 

• The LSP momentum cannot be reconstructed except for a few very special part of 
the decay phase space. 

• Only the events near the endpoints of the decay distribution are used for the mass 
fit. The large statistics are required to see the endpoints while the events away from 
the endpoints also contain the independent information of the masses. 

• The selected SUSY events may contain events from several cascade chains. For 
example, three or four channels involving hi or tj contribute simultaneously for hhll 
final state. The measured m{hhll), m{hll) and m{hl) endpoints are the weighted 
average of the endpoints of these channels. This introduces additional systematical 
uncertainty for the mass determinations. 



tan (3 


m-g 










10 


595.2 


491.9 (524.6) 


176.8 


143.0 


96.0 


15 


595.2 


485.3 (526.9) 


177.9 


143.0 


96.5 


20 


595.2 


478.7 (531.2) 


178.5 


143.1 


96.7 



Table 2: Masses of some sparticles at SPSla and the variants studied in Section 3, in GeV. 

The above problems may be removed by using the "mass relation method" |12|13j 
when a sparticle cascade decay chain is sufficiently long. In the mass relation method, 
one uses the on-shell condition of the sparticle masses to solve the kinematics of the 
sparticle decay product exactly, and reconstruct the masses of the SUSY particles as the 
peaks of certain distributions. To describe the method, let us take the process 

-g^hh^ X%h ^ ihh£2 ^ XibMh, (7) 

where 6i, 62, ^1 and I2 denote different b quarks and leptons respectively. Both of the 
sbottom states bi and 62 yield the decay chain of Eq. ((Tj). 

Five sparticles are involved in the cascade decay Eq. (|7j), therefore one can write five 
mass shell conditions among the leptons and quarks in the final decay products. 

^? = 

= (hi + + +PHf, 

= (Px'i + P^i + P^2 + Pbi + Pbif ■ (8) 

For a bbii event, the equations contain the 4 unknown degrees of freedom of the 
Xi momentum. Each event therefore describes a 4-dimensional hyper-surface in a 5- 
dimensional mass parameter space, and the hyper-surface differs event by event. From 
the purely mathematical point of view 5 events would be enough to determine a discrete 
set of solutions for the masses of the involved sparticles. 

From now on, we will develop the argument by assuming that the masses of X2? ^) 
and Xi a^^e known for simplicity. This is a reasonable assumption at the LHC, where 
it has been shown that a detailed study of the lepton-lepton system from the X2 decay 
can be used to precisely constrain these masses. In addition to that, if the international 
linear collider (ILC) would be build, the mass of lighter sparticles will be measured with 
less than 1% accuracy there. If we assume the mass of the lighter sparticles is known, 
each event corresponds to a different curve in the (m^, m^) plane, which is expressed as 
following equation; 

f {nig, mi, pbi, Pb2, Pii, P12) = Qiim'l + 2Qi2mlml + Q22mf + 2Qiml + 2Q2ml + Q = 0, (9) 

where Q's are the functions of the lighter sparticle masses and momenta of the b quarks 
and leptons. Two events are enough to solve the gluino and sbottom masses altogether 
up to the maximally four solutions. ^ 



^In addition, there are two possibihty to assign the two leptons to li and l2- In the plot of this paper, 
we fix the assignment by fixing the assignment so that the lepton with the higher px is h. 
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Figure 1: The gluino mass distributions 
for SPSla with the event pair analysis. 
The open, green, and blue histograms 
in the top figures are for OSSFxOSSF, 
OSSFxOSOF, and OSOFxOSOF event 
pairs, respectively. The open histograms 
in the bottom figures shows the mass dis- 
tributions after background subtraction. 
The contributions of 62 are shown by red 
histograms. 



We call this technique "mass relation method" , 
because here one uses the fact that the sparticle 
masses are common for all events which go though 
the same cascade decay chain. Note that the events 
need not to be close to the endpoint of the decay 
distribution, but they are still relevant to the mass 
determination. This means that one can use the 
mass relation method even if the number of signal 
events is small. 

In the following, we start from the selected bbii 
events at SPSla. The decay distribution is stud- 
ied by generating events by using HERWIG^^, and 
simulate the events by using the ATLFAST detector 
simulator In order to minimize the combinato- 
rial backgrounds we use the event pairings which 
satisfy the following conditions. 

• Eq. has solution for only one of the two 
possible lepton assignments. 

• For the selected lepton assignment the result- 
ing quartic equation in m? has only two solu- 
tions, and the difference of the gluino masses 
for the two solutions is more than 100 GeV. 
The smaller gluino mass solution is chosen. 

Note the selections are rather phenomenological and 
they may introduce some bias to the reconstructed 
sparticle masses. 

The rrig distributions for the OSSFxOSSF 



events pairs are shown in the histograms on the up- 
per line of Fig. [T] We find the obtained peak po- 
sitions are consistent with the input masses in this 
study. A significant SUSY background, also shown in Fig. ^ is still present in the sample. 
This background can be estimated from the data themselves by using the bbii events with 
an opposite sign opposite flavor (OSOF) lepton pair (i.e. ii = e^ix^). 

The histogram shows a peak corresponding to the input value for the gluino mass even 
before the background subtraction. The green and blue histograms show the estimated 
background distributions. The distributions after the background subtraction are shown 
in the histogram on the lower line of Fig. ^ The peak position and its error obtained by a 
Gaussian fit to the distribution are and it is 591.9±0.7 GeV at SPSla. One can also look 
into the distribution of — and estimating the value of this observable by performing 
a Gaussian fit on the observed peak and the peak position is 98.9 GeV, while the input 
is 103.3 GeV. 



The error of the peak position is 0(1) GeV, but this is not the error of the gluino 
mass, because each event can be used many times for the mass calculation. To estimate 
the error, we generate many signal events at SPSla(tan/3 = 10) and study the deviation 
of the peak positions of the sub-samples. The gluino mass error for / dtC = 300fb~^ is 
Anig ~ 2 GeV at SPSla. 
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Figure 2: The distributions for tan/3 = 10 (left), 15 (center) and 20 (right) with a fixed gluino mass 
{rrig = 595 GeV). The open and green histograms in the top figures show the distributions of OSSF and 
OSOF lepton events, respectively. The mass distributions after background subtraction are shown in the 
bottom figures, where the contribution of the 62 events are shown by the red regions. 

Once the gluino mass is fixed by the analysis shown above, Eq. (0) can be solved for 
each event for two solutions for the sbottom mass, giving as input the central measured 
value for the gluino mass. In Fig. we plot the distribution of the smaller sbottom mass 
solution m^(min) for both OSSF (signal) and OSOF (background) lepton pair events (top 
histograms) for SPSla and two other model points which has same MSUGRA parameters 
but tan/5 = 15 and 20. The relevant sparticles masses are listed in Table □ The m^^ 
changes sensitively with tan/3 due to the increased left right mixing. We show in the 
bottom line the mass distributions after the background subtraction. The peak positions 
are evaluated by a Gaussian fit, and listed in Table El Note that the total number of the 
signal bi event is smaller by factor of 4 for tan/? = 20 compared with tan/? = 10. This is 
because the decay X2 ~^ '''^ dominates as tan j3 increases. Thanks to the peak structure of 
the distribution, the mass peak is still seen very clearly. The mj^(min) peak and m^(input) 
are in good agreement, and it is not so for the larger solution m^(max). 

The peak positions of the distribution of the events originated from 62 decay are also 
consistent with 62 masses. The existence of 62 can be established only after understanding 
b jet smearing and bi distribution correctly. 
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Table 3: Fit results of the sbottom mass in GeV with a fixed gluino mass {nig — 595 GeV). 



2.4- Likelihood analysis 

Because Eq. Q have multiple solutions for a fixed gluino mass, we encountered the 
problem to select the one of the multiple solution artificially in the previous subsection. 
To avoid bias to the analysis, we introduce a "likelihood analysis" in this subsection. 

From Eq. 0, each event is represented as a curve in the (mg,m^) plane. The coef- 
ficients of the curve are a function of the four momenta of the detected partons. The 
partons are measured as jets in the detector, which smears the parton according to a 
smearing function. From the measured quadratic form for an event we can build a con- 
fidence belt in the {mg,mi) plane [TH|13j . In order to build the "probability distribution 
function", the crucial ingredient is the distribution of the measured 6-jet momenta as a 
function of the 6-parton momenta. We define an approximate probability density function 
according to the formula: 

C(m„m,) =l,p' jdp'AP. p'Mp. ■ pWnm„,nM) (10) 

where e{p : p') is the probability to measure a momentum p' for a b jet, given a 6-parton 
with momentum p. 

In Eq. fll()|l we did not include the possibility of lepton momentum mis-measurement, 
which has an almost negligible effect. We also assume that the jet direction is not modified 
by the measurement and we use for e{p : p') a gaussian distribution, with a width a 
corresponding to the parameterized jet smearing used in the fast simulation program. The 
approximate function takes however into account the dominant part of the jet smearing 
and can be used to demonstrate the method. 

We now show numerically calculated log C in the {rrig — mi, nig) plane for a few events 
where the hhCC events originates from the cascade decay of Eq. ((Zj) at SPSla. We calculate 
£ using the following procedure. For each event in our sample, characterized by a (^1,^2) 
pair of measured momenta for the 6-jets, we generate Monte Carlo events of two h jets with 
momentum (p'^, ^2), where p'^ and p'2 are randomly generated according to the function 
e(pi : p'l) X e{p2 : ^2)- The histgrams of the number of curves which satisfy Eq. that 
go through 1 GeVxl GeV grid in the [nig, mi) plan, normalized by dividing the total 
generated event n corresponds to C{mg,mi). In Fig. 01 we plot 

Alog£ = \og{C{mg,mg + Amg,m^,m^ + Am^) + c) - log(£(min)), (11) 

where c = 0.001 is a constant cutoff factor, which is needed as for each event we generate 
only a finite number of Monte Carlo experiments, and therefore some bins can have 




zero hits. The shape of the probabihty density distribution is different event by event, 
as it depends on the event kinematics. The size of the allowed band is also different, 
which means that some events will have more weight in the determination of the mass 
parameters. 

By combining the probabilities for different events, a region of maximum probability 
in the {nig, mi) is found, where the curves of maximum probability for all events approx- 
imately cross. Namely, we can build the combined likelihood for all the events defined 
as: 

log£comb("^g,"^g + Amg,m^,m^ + Am^) = ^ log(/:(mg, + Am^, m^, + Am^) + c). 

events 

(12) 

To study the likelihood distribution of gluino and sbottom mass, we actually shows the 
distribution of subtracted likelihood 

log£sub = log^ossF -log^osoF = ^ogC- ^ log£. (13) 

OSSF OSOF 

This is not the correct definition of the likelihood function, but in the limit of infinite 
statistics, log£sub is independent from the contribution of accidental lepton pairs. For 
the correct treatment, see [T^ . 

We plot the contours of the function log £sub in Fig. IH where plots (a) and (b) [(c) and 
(d)] are for tan/3 = 10 [tan/3 = 20]. The distributions (a) and (c) are produced accepting 
all the events which pass the selections, whereas distributions (b) and (d) are produced 
using an event sample where the events including a 62 decay have been rejected. 

In Fig. 0] (a) and (c), the position of the peak for irig — is roughly consistent with 
the input value. Unlike the gluino and sbottom mass fits in the previous section, we 
obtain the correct peak position without the need of artificially choosing among multiple 
solutions. The likelihood distribution can be used to determine the g and h. We restrict 
the likelihood distribution for 591 GeV < uig < 599 GeV (within 4 GeV from the input 
gluino mass). We then fit the distribution around the peak assuming gaussian distribution. 
The likelihood distribution peaks at the gluino and sbottom mass difference as 99.5 GeV 
for tan/3 = 10, 104.2 GeV for tan/3 = 15, and 113.9 GeV for tan/3 = 20, where the input 



(a) tan|S=10, all (b) tan|3=10, no sbottom2 (c) tan/S=20, all (d) tan|S=20, no sbottom2 
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Figure 4: Contours of the likelihood function log£sub in the (mg — mi, rrig) plane: (a) for (c) tan/3 = 10 
and (b) and (d) for tan/3 = 20, respectively. The contours (b) and (d) are made without 62 contributions. 
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Figure 5: The likelihood as a function of rrig — for tan/3 = 20: (a) for all events and (b) without 62 
events. 

value is 103.3 GeV, 109.9 GeV and 116.5 GeV, respectively. The fitted values display shift 
of about 4 GeV from the true value. We ascribe this effect to our simplified modeling 
of the jet smearing in building the likelihood function, which should disappear once the 
detector response is properly taken into account in the unfolding procedure. 

In Fig.El^a), we show the distribution of log £sub as a function of rrig— mi at tan [5 = 20, 
restricting the gluino mass in the region 591 GeV< rrig <599 GeV again. On the left of 
the peak corresponding to the 61 mass, we see a small bump in the distribution. This 
bump is not observed in the mass distribution made without 62 contribution (Fig. |31|^b)). 
In order to claim the presence of a second component in the distribution on the data, the 
ability of correctly reproducing the likelihood distribution for bi events would be needed. 
It is also difficult to extract a statistical significance for the 62 shoulder as our definition 
of the likelihood function is approximate one. 

2.5. More challenges 

In the previous subsection, we have discussed the MSSM parameter determination 
at SPSla. For this point, the mass spectrums would be obtained quite accurately by 
looking into the decay cascades involving leptons. The SUSY scale is relatively light 

~ '^q ~ 600 GeV, when production cross section is as large as 50 pb. One can even 



determine the edges of the x1 decays. If the squark and gluino mass scale is around 1 
TeV, the production cross section is only around 3 pb, and statistical errors would be 
increased accordingly. 

The sensitivity to the masses would be further reduced when rriq > rrig. For this case 
squarks decay dominantly into gluino, and the gluino decay is 3 body. Then it is hard to 
identify the jets from squark and gluino decays. This is typically the case for the model 
points in the "Focus point region" in the supergravity model. 

The mass resolution also may be worse when tan (3 is large in MSUGRA, because the 
branching ratios of the X2 i^^o n or e are reduced significantly as Br{x2 ~^ ^'^) dominates. 
The decay into r may also be studied, but a r jet resolution is significantly worse compared 
to that of a lepton, and the background are much higher due to fake r's from QCD jets. 
The case where f and Xi are degenerated in mass is favored in MSUGRA model with the 
dark matter constraint, because the coannihilation of f and Xi reduces the thermal relic 
density in the Universe to the acceptable level. In that case, one of the r lepton from X2 
decays may be too soft to be detected at the LHC. 

It is therefore important to establish analysis which does not relay on the leptons in 
the event cascades. For this purpose we discussed the reconstruction of the tb final state 
where top decays hadronically [T7I . The dominant decay of gluinos for these points are 
g {i,b)-^ tbxf- The top quarks in the SUSY events may be reconstructed by looking 
for the bjj where rrijj ~ myi/ and m^jj ~ m^. There are many accidental background 
which satisfy this conditions, because the number of jets njets in the SUSY events are 
typically 8 to 10. Such backgrounds can be subtructed by using the side band events 
where nijj < m^y — 15 GeV and rrijj > mw + 15 GeV. The reconstructed top quarks are 
then used to study the nitb distributions. In Fig. El we show the rritb distribution after 
the background subtruction for SPSla and SPS2, where SPS2 corresponds to the focus 
point where rriq ~ 800 GeV and ~ 1.5 TeV. For SPSla, g tt and bb is open, and 
distribution shows sharp edge consistent to those calculated from the input masses. The 
endpoint is measured within an error Am^f, ~ 4 GeV, and the distribution maybe used 
to study the t nature. For the case of SPS2, the decay chain is not open, but gluino still 
dominantly decay into third generation particles because stop and sbottom masses are 
significantly lighter than that of the first and second generation squarks. The distribution 
does not show the edge strucuter because the gluino decay is three body. The rritb endpoint 
is consistent to the mass difference — m-+ ~ 560 GeV and — m-+ ~ 480 GeV. 

3. Supersymmetric Relations 

In MSSM, all dimensionless couplings involving supersymmetric particles are described 
by either SM Yukawa or gauge couplings. Not only the the size of the couplings, but 
the chiral structure of the couplings are restricted by the supersymmetry. The fermion- 
sfermion-gaugino coupling is restricted to preseve the chirality of (s)fermion, therefore for 
example. 



WlJl + tan OwBIlII - 2 tan^vK^/i?/^ + h.c. 



(14) 



V2 

for the lepton-slepton-gaugino interaction, while the chirality must be flipped for the 




Figure 6: The mtb distribution at SPSla and SPS2. The formar is dominated by the two body cascade 
decay g ti, bb tbx^ and the latter is dominated by the three body decay g tbxf ■ 

Higgsino interactions. Proving the couplings and their chiral structures is important 
steps to check if the discovered new particles are indeed superpartners of the SM particles 
or not. 

For some cases one can check the chiral structure of the interaction by looking into 
the decay distributions. We first note that the SUSY particles are naturally polarized in 
many SUSY process. For example, in the cascade decay 

Ql xW IrIrQ llqXi, (15) 

the X2 is polarized as right-handed (opposite to ql), because the Yukawa coupling of 
squark- quark- ino flip the chirality. The polarized X2 further decay into either leptons- 
antislepton, or antilepton-slepton. The two decay branching ratios are same because X2 
is a Majorana particles. For the decay X2 ^ ^bJ-~^ through Yukawa type interaction, anti- 
lepton likely to go in the same direction to X2 (namely opposite to the jets in the squark 
rest flame), while the lepton goes in the same direction to the jet for the X2 ~^ ^*r^~ decay. 

The difference of the angular distribution appears as the charge asymmetry in the 
m{jl) distribution |18j. becuase the m{jl) is propotional to 1 — cos 6', where 6 is the angle 
between the jets and the lepton in the X2 ^^st flame. The m{jl^) is peaked sharply 
at the endpoint of m{jl) distirubtion, while the m{jl~) distribution is suppresed at the 
endpoint. Finally a (anti) lepton is emitted from a (anti)slepton, but this time the lepton 
angle distriution is spherical in the I rest flame. 

The charge asymmetry is exactly opposite for q1 decays. However, at pp colliders, 
the number of produced squark is much larger than that of anti-squark, and the X2 
produced from qi decay dominates total X2 productions. The charge asymmetry in m{jl) 
distribution in the jU sample therefore remains. From the distribution, one can conclude 
that the decay products X2 is a fermionic particle, and the interaction is chiral [19^20] . 

The detectability of the charge asymmetry can be studied by generating events by us- 
ing HERWIGjl4j. and simulate the events by using the ATLFAST detector simulator |15j. 
Here we show two representative case at the points SPSla and SPSS in Fig. [71 where 
the decay X2 ^ ^^R dominates at SPSla and Ul dominates at SPSS respectively. The 
distributions shows opposite asymmetry near the endpoints the m{jl) distributions for 



the decay q ^ qX2 ~^ qll (298 GeV at SPSla and 200 GeV at SPSS respectively). This 
is because the chirahty of the X2"^"^ couphng is exactly opposite for the two points [ 



SPS1a(tan|3=10) 



200 400 
m(ll)[G8V] 



Figure 7: The reconstructed asymmetry [Nsig{il+) - iVsig(j7")]/[7Vsig(j7+) + 7Vsig(j7")]. 

The cascade decay may also be sensitive to the II and Ir mixing PU]. Due to the F 
term contribution of slepton mass matrix, ~ m//i tan/?, the lighter slepton mass 

eigenstates li is mixture of II and Ir where the mixing angle is defined as li = II cos 6'; + 
iRsmOi. The left and right hand coupling of the to the lighter mass eigenstate li is 
expressed as 

C = -^xI{ClPlI + CrPrI)11, 
Cl = -cosei{N^^ + N*^^t8inew), 

Cr = 2tan^iy sin^/A^^2> (16) 

when X2 is gaugino like. 

When X2 ~ Xi ~ -S, and h ~ Ir, Cl is suppressed by small mixing angle pro- 
portional to rriffi tan (3, while Cr is suppressed by the small bino component of the wino 
like X2- The non- universality appears as the sign of the non-zero F term mixing of ft and 
visible in the wide parameter region. 

While the mixing angle of fi is only around 1% for the Snowmass point SPSla, the 
difference of the decay branching ratio Br{x2 /^A'')/-Br(x2 c/^) from 1 is as large as 
4%. Although the statistics reduces as tan/? increases, statistical significance of the signal 
defined as 

/fir(^_ \ / AiV(sig) \-^ 
[Br{e) ) \ iV(sig) ) 

where AA'"(sig)= \/Nos, and A^^os is the number of the total odd sign two lepton events. 
We found that S increase from 5.4 at SPSla to 8.5 for a modified point with tan/5 = 20 
(see Table El). 

4. Discussion 



In this talk, we discussed SUSY studies at the LHC In early '90, the LHC was con- 
sidered merely as a discovery machine, due to the challenging experimental environment. 



point 


Br{x2 ee) 


iVsig(e and /i) 


A^os 


\Br(e) ^) 


s 


tan /3=10 


6.3% 


1.39 X 10^ 


2.68 X 10^ 


4% 


5.4 


tan (3 = 20 


1.2% 


0.28 X 10^ 


1.02 X 10^ 


17% 


8.5 



Table 4: The Br{x2 — > ee) and accepted number of events for SPSla tan (3 = 10, 20. iVgig is the number 
of events after e^/i^ subtraction, while A^os is the number of the total odd sign two lepton events. 
The expected deviation of Br{x2 ~* Mm) from that for ee is compared with the statistical error of the 
number of the e+e~ events for J dtC = 300 fb~^. 

However, it was recognized that the parameters of the MSSM can be determined at the 
LHC Since then, many techniques have been developed to extract the mass and couphng 
information. 

When the gravitino is the LSP and the NLSP is long lived, the LHC works as the 
machine for precision studies. This is because SUSY events contains many leptons and 
photons whose momentum can be easily combined to solve sparticle masses. For the case 
of the charged NLSP, the semi-stable charged particles would be observed as a highly 
ionizing charged tracks which goes through the detector. The decay position and decay 
time of the NLSP can be measured at the LHC, and the gravitino interaction to the 
matter may be explored. The cosmology related to the gravitino LSP and the collider 
phenomenology have received lots of attentions. 

When the LSP is the lightest neutralino, the sparticle mass determination using the 
endpoint analysis is known to be very powerful. For the most favorable case, one can 
determine the LSP mass within 5% and squark and gluino masses within 2%. We have 
also introduced the mass relation method developed very recently in this talk. By using 
this method, one can reconstruct sparticle masses as peaks in the distributions calculated 
from the events using the exact formula to solve decay kinematics. The mass relation 
method may be useful even if the event statistics is small. 

Finally we discussed the study of the sparticle interaction at the LHC. The chiral 
structure of the sparticle couplings may be studied through the charge asymmetry of the 
selected m{jl) distributions, where the jet and the lepton come from the cascade decay 

ql xlq llq- 

To complete the task to explore all masses and couplings of SUSY particles is probably 
impossible by the LHC alone. However even after the ILC is started, the measurements 
at the LHC are useful to understand supersymmetry in nature. There are many on-going 
studies, and we still need more ideas and thoughts so that we do not miss any new physics 
signatures at the LHC. 
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